Fast Equilibration of Hadrons in an Expanding Fireball 
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Due to long chemical equilibration times within standard hadronic reactions during the hadron 
gas phase in relativistic heavy ion collisions it has been suggested that the hadrons are "born" 
into equilibrium after the quark gluon plasma phase. Here we develop a dynamical scheme in 
which possible Hagedorn states contribute to fast chemical equilibration times of baryon anti-baryon 
pairs (as well as kaon anti-kaon pairs) inside a hadron gas and just below the critical temperature. 
Within this scheme, we use master equations and derive various analytical estimates for the chemical 
equilibration times. Applying a Bjorken picture to the expanding fireball, the kaons and baryons 
as well as the bath of pions and Hagedorn resonances can indeed quickly chemically equilibrate for 
both an initial overpopulation or underpopulation of Hagedorn resonances. Moreover, a comparison 
of our results to (B + B)/ir + and K/ir ratios at RHIC, indeed, shows a close match. 

PACS numbers: 



(Anti-)strangeness enhancement was first observed at 
CERN-SPS energies by comparing anti-hyperons, multi- 
strange baryons, and kaons to pp-data. It was consid- 
ered a signature for quark gluon plasma (QGP) because, 
using binary strangeness production and exchange reac- 
tions, chemical equilibrium could not be reached within 
the hadron gas phase 0|. It was then proposed that 
there exists a strong hint for QGP at SPS because 
strange quarks can be produced more abundantly by 
gluon fusion, which would account for strangeness en- 
hancement following hadronization and rescattering of 
strange quarks. Later, multi-mesonic reactions were used 
to explain secondary production of anti-protons and anti- 
hyperons 0,0. At SPS they give a typical chemical 
equilibration time Ty « 1 — 3^r using an annihilation 
cross section of a pY ~ o~ p p w 50mb and a baryon density 
of pb ~ Po to 2po, which is typical for SPS. Therefore, 
the time scale is short enough to account for chemical 
equilibration within a cooling hadronic fireball at SPS. 

A problem arises when the same multi-mesonic reac- 
tions were employed in the hadron gas phase at RHIC 
temperatures where experiments show that the parti- 
cle abundances reach chemical equilibration close to the 
phase transition 0]. At RHIC at T = 170 MeV, where 
a ~ 30mb and p C g w p e g w 0.04fm~ 3 , the equilibrium 
rate for (anti-)baryon production is r « 10—, which 
is considerably longer than the fireball's lifetime in the 
hadronic stage of r < 5=jp. Moreover, r as 10^p was also 
obtained in Ref. H using a fluctuation-dissipation theo- 
rem and a significant deviation was found in the popula- 
tion number of various (anti-)baryons from experimental 
data in the 5% most central Au-Au collisions [fjj. These 
discrepancies suggest that hadrons are "born" into equi- 
librium, i.e., the system is already in a chemically frozen 
out state at the end of the phase transition [7J, |8| . 

In order to circumvent such long time scales it was 
suggested that near T c there exists an extra large par- 



ticle density overpopulated with pions and kaons, which 
drive the baryons/anti-baryons into equilibrium [9(. But 
it is not clear how this overpopulation should appear, and 
how the subsequent population of (anti-)baryons would 
follow. Moreover, the overpopulated (anti-)baryons do 
not later disappear [l(|. Therefore, it was conjectured 
that Hagedorn resonances (heavy resonances near T c with 
an exponential ma ss sp ectrum) could account for the ex- 
tra (anti-)baryons [lfj| . Baryon anti-baryon 1^, U| and 
kaon anti-kaon production develop according to 



mr <-> HS <-» rii^TT + BB, 
rn: <-* HS «-> n,i,fc7r + KK, 



(1) 



which provide an efficient method for producing baryons 
and kaons because of the large decay widths of the Hage- 
dorn states. In Eq. ([T]), n is the number of pions for 
the decay nir <-> HS, which can vary, and (ii,fe) is 
the number of pions that a Hagedorn state will decay 
into when a baryon anti-baryon (kaon anti-kaon) pair is 
present. Since Hagedorn resonances are highly unstable, 
the phase space for multi-particle decays drastically in- 
creases when the mass is increased. Therefore, the reso- 
nances catalyze rapid equilibration of baryons and kaons 
near T c where the Hagedorn states show up. Here we use 
a Bjorken expansion within a cooling fireball in order to 
see at which temperature the chemical equilibrium values 
are reached or maintained. In this letter we also briefly 
discuss an analytical solution of the chemical equilibra- 
tion time, which is valid at a constant temperature near 
T c . Moreover, our numerical results for the baryon anti- 
baryon pairs and kaon anti-kaon pairs suggest that the 
hadrons can, indeed, be born out of equilibrium. 
We use a truncated Hagedorn mass spectrum [l3j 
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where the Hagedorn temperature is set to Tu = 180MeV, 
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which lies within the present range of Lattice QCD pre- 
dictions 14 1, the normalization factor is A — 0.5 MeV 5 , 
and mo = 0.5 GeV. We consider only mesonic, non- 
strange resonances and discretize the spectrum into mass 
bins of 100 MeV that range from the mass M = 2 GeV 
to M = 7 GeV. The effects of the truncation and Hage- 
dorn temperature are further discussed in [T^. However, 
the values we have chosen are acceptable. 

The abundances' evolution of the Hagedorn states, pi- 
ons, and baryon anti-baryon pairs due to the reactions in 
Eq. (fTJ) are described by the following rate equations 



where J ^dy = 874 from Ref. [H for the 5% most cen- 
tral collisions within one unit of rapidity and the entropy 
per pion S^/N^ = 5.5 is larger than that for a gas of 
massless pions according to the analysis in Ref. [3] . The 
volume [3j is 
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(5) 

where the initial radius is ro(to) = 7.1 fm, the av- 
erage transversal velocity varies vq — 0.3c, 0.5c, and 
0.7c, and the corresponding accelerations are taken as 
a = 0.035,0.025, and 0.015, respectively, 
j The equilibrium values of pions, N% q , shown in Fig. [1] 
are found using a statistical model |20j . Here we consider 
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where the fugacity is A = p, N is the total number 
of each particle, and its respective equilibrium value is 
N eq . The summation over i represents the i th Hagedorn 
resonance bin. The structure of the rate equations for the 
kaon anti-kaon pairs is the same as in Eq. (J3]), however, 
K K is substituted in for BB. 

The branching ratios for HS <-> mr are described by 

("-<"i>) 2 

a Gaussian distribution Bi „ « — i=e 2cTi where 
(rii) = 0.9 + 1-2^- is the average pion number each 
Hagedorn state decays into, found in a microcanonical 
model (l5j]. of = (0.5^- ) 2 is the chosen width of the dis- 
tribution, and n > 2 is the cutoff for the pion number. 
Moreover, the branching ratios are normalized such that 
Y^=2 Bi,n = 1, which gives (rij) » 2 to 9 and of w 0.8 to 
11. The total decay width, r f « 0.17mj-88 MeV, which 
ranges from r.; = 250 to 1100 MeV, was found using the 
mass and decay widths in [l6j and fitting them linearly 
similarly to what was done in Ref. [l7| • The decay widths 
for the baryon anti-baryon decay are T i BB = {B)Ti and 
Ti i7V = Ti — T { BB . The average baryon number (B) per 
unit decay of Hagedorn resonances within a microcanon- 
ical model ranges from 0.06 to 0.4, so T i BB — 15 to 400 
MeV [UJ. We use only the average values in Eq. fl} so 
that (rii t b) = (ni,k) = 3 to 6 1(| 15 1 is used for both the 
baryons and kaons. For the kaons K ^ = (K)Ti where 
(K) = 0.4 to 0.5 [HI [IB)]. Thus, heavier resonances equi- 
librate more quickly because of large decay widths. 

Using a Bjorken expansion, we find a relationship be- 
tween the temperature and the time, T(t), for which the 
total entropy is held constant 
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FIG. 1: Comparison of the total equilibrium number of pi- 
ons iVp?, Hagedorn states Nf 9 , and effective pions N% q as 
defined in Eq. ©. 

both the direct pions and the indirect pions, which come 
from resonances such as p, u> etc, and both the direct and 
indirect kaons. In Fig. Q] we see that N^ q increases with 
decreasing temperature. This occurs because the Hage- 
dorn states dominate the entropy at high temperatures, 
which affects N^ q due to the entropy constraint in Eq. 
(f?|). Therefore, we must consider the number of "effec- 
tive pions" in the system, i.e., the total number of pions 
plus the potential number of pions from the Hagedorn 
resonances, defined as 



-jhim) + — p — (n»,fc) 



K,bb = ^ + ^V i f^(r, i } + ^K b )j (6) 

for the kaons and baryons, respectively. In both cases 
N% q remain roughly constant throughout the Bjorken ex- 
pansion. Additionally, throughout this paper our ini- 
tial conditions are the various fugacities a = X v (to), 



Pi 



Xi(t ), and 4> = X BB (t ) or = X KK (t ), 



which are chosen by holding the contribution to the 
total entropy from the Hagedorn states and pions 
constant i.e. s Had (T , a) V(t ) + sjjg(T , Pi)V(t ) = 
SHad+Hs(To)V(t ) = const. 
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FIG. 2: Numerical and analytical results in a static environ- 
ment at T — 180 MeV when # = 1.3 and a = 0.7. 
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FIG. 3: Same as Fig. [2] with no initial baryons (<f> = 0). 



The initial estimate for the Hagedorn state equilibra- 
tion time is r,; = 1/I\. In order to estimate the chemical 
equilibration time, we use Eq. (3} in a static environ- 
ment to find the equilibration time to be in the general 
ballpark @, [T3] of 
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between T = 180 to 170 MeV. As will be proven in detail 
in Ref. 12j, these time scales are only precise when the 
pions and Hagedorn states are held in equilibrium. In 
reality the chemical equilibration times are more compli- 
cated due to non-linear effects and the evolution of the 
equilibration must be divided into separate stages for a 
sufficient analysis. 

To find time scale estimates, we consider the more sim- 
plified case near T c excluding the baryons and kaons, i.e., 
Eq. without the baryonic terms. The evolution of the 
rate equations can be divided into three stages as shown 
in Tab. Q] and derived in [13]. Initially, when the pions 
are far from equilibrium (X n rs 0) the Hagedorn states 
can be held constant at a constant fugacity Substi- 
tuting Att ~ and A-j ~ (3i into Eq. (J3]), we obtain t°. 
As the pions near equilibrium, we can then use A^ — > 1 
to obtain t„. Eventually, the right-hand sides of Eq. ([3]) 
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TABLE I: Equilibration times from analytical estimates 
where QE is quasi-equilibrium and TOT is total equilibrium 
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FIG. 4: Results for an expanding fireball when a = 1, f3 = 1, 
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and <f> = 0. The effective number of baryons kaons 
NVtf,, and pions N eq ai , and N eq 

KK 1 c 7T.BB 



become roughly zero before full equilibration (known as 
quasi-equilibrium), which occurs once the lightest reso- 
nance reaches quasi-equilibrium v 2GeV — 0.8^p. To ob- 
tain t® e we solved Eq. © without the baryonic term, 
assuming A x — > 1 and that the right-hand side of the pion 
rate equation equals zero. Then the total equilibration 
(7) time r tot is just the addition of TiGeV and t® e . Since r tot 
includes all the non-linear effects, which occur even after 
equilibrium is neared, the more appropriate time scale 
is on the order of t®. The equilibration times increase 
directly with N% q , (nf) and are shortened by large IYs 
and wide branching ratio distributions Cj's. Because N^ q 
decreases quickly as the system is cooled, the equilibra- 
tion time is significantly longer at lower temperatures. In 
Fig. [2] our analytical fit, which are exponential fits [12] 
based on Tab. [H match our numerical results well and 
nicely concur with the numerical results in Fig. [3] Addi- 
tionally, the baryons take slightly longer than predicted 
in Eq. ([7]), but they still equilibrate quickly (Fig. [3]). 

In Fig. [4] the baryons and kaons are shown for an ex- 
panding system where we see that the baryons reach 
chemical equilibrium by T = 165 MeV (t-t «2-3^) 
and the kaons at T = 160 - 140 MeV. As with the pions, 
we consider the effective number of baryons and kaons 
because of the effects of Hagedorn resonance on the en- 
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FIG. 5: Comparison of the ratios obtained employing the 
expanding fireball picture for various initial conditions to data 
at y^s = 200 GeV from PHENIX [U and STAR [H . 



tropy at high temperatures, so 



N BS = N 



BB 
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which are shown in Fig. [4] Not surprisingly, N^ q gg and 
N eq T , r> remain almost constant due to the constraint set 
in Eq. Moreover, our expansion is not strongly af- 
fected by Vq and, therefore, in the following graphs it is 
set to vq = 0.5c. 

In Fig. [3] we compare our total baryon to pion ratio 



21] 



to experimental data from PHENIX 
2j. (B + B)/tt+ is calculated by B + B = 
n + n « 2(p + p). It should be noted here that 



(B + B)/ir+ 
and STAR " 
P + P 

in our calculations we use both the effective number of 
baryons, in Eq. ([5]), and pions, in Eq. ©. We obtain 
(B + B)/it + w 0.3, which matches the experimental data 
well. Moreover, our results are independent of the chosen 
initial conditions. Also, in Fig. [5] we compared the kaon 
to pion ratio to the data at PHENIX [2lj] and STAR [22 1 
(both K/tt + and K/tt + are shown). Again, we use the 
effective number of kaons (J8]) and pions ([6]) . Our K/n + 
ratios compare to the experimental data very well and 
they level off between 0.16 to 0.17. As with the baryon 
anti-baryon pairs we do not see a very strong dependence 
on our initial conditions. In Fig.[5]both figures agree well 
with experimental data. Moreover, they remain roughly 



constant after T = 170 — 160 MeV. This demonstrates 
that the potential Hagedorn states can be used to explain 
dynamically the build up of the known particle yields. 

In future work we will consider strange baryonic de- 
grees of freedom and thoroughly study the effects of 
our initial conditions and parameters. To conclude, we 
used Hagedorn resonances as a dynamical mechanism to 
quickly drive baryons and kaons into equilibrium between 
temperatures of T = 165 — 140 MeV. Once a Bjorkcn 
expansion was employed, we found that our calculated 
K/n + and (B + B)/ir + ratios matched experimental data 
well, which suggests that hadrons do not at all need to 
start in equilibrium at the onset of the hadron gas phase. 
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